The nature of an enhanced ferroelectric phase transition temperature in 

perovskite-based solid solutions 
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We explain the phenomena of ferroelectric phase transition temperature T c enhancement beyond 
the end members in perovskite solid solutions like BiMeC>3-PbTi03 (Me=Sc, In, etc.) is related 
to nonlinear and spatial correlation effects. The explanation is based on the calculation of T c 
in the framework of our random field theory with additional account for nonlinear effects in the 
above substances. We show that the maximum of T c for certain PbTiOa content takes place when 
coefficient of nonlinearity is positive, the value of this coefficient is found from best fit between 
theory and experiment. This nonlinearity coefficient is the only adjustable parameter of the theory. 
We show that enhancement of positive nonlinearity coefficients enhances greatly the T c maximum 
over its value for end members. 

The theory lays the foundation to calculate not only T c for above solid solutions but virtually 
any equilibrium and/or nonequilibrium thermodynamic characteristics such as static and dynamic 
dielectric susceptibility, specific heat etc as a functions of PbTiC>3 content, temperature, electric 
field and other external parameters. 

PACS numbers: 77.80.Bh,77.84.Lf,77.84.Dy 
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I. INTRODUCTION 

Over the last few years considerable effort has been 
spared to synthesize the dielectric materials with con- 
trollable properties for many technical applications. The 
most promising substances are believed to be the com- 
pound materials consisting of solid solutions of differ- 
ent combinations of ferroelectrics with different dielec- 
tric properties. One of the examples is perovskite solid 
solutions like BiMe0 3 -PbTi0 3 (Me=Sc, In, etc.), which 
have high ferroelectric phase transition temperatures T c 
at the morphotropic phase boundaries with enhancement 
beyond the end members. Such materials can be used as 
the materials with excellent high-temperature piezoelec- 
tric propertie»ii2ii^i&. 

A common feature of the above solid solutions is the 
existence of numerous random fields sources due to sub- 
stitutional disorder, unavoidable impurities, vacancies in 
anion and cation sublattices etc. These random fields 
play a crucial role in the properties of disordered ferro- 
electric and magnetic materials (see e.giS*^) and above 
substances in particular. This means that observable 
physical properties of above systems depend strongly on 
the form of random fields distribution function. Namely, 
the relation between width of the distribution function 
(the dispersion of random fields) and its first moment 
(mean value of random field) generates all observable 
features of the phase diagram of a disordered dielectric 
and/or magnetic materials, i.e. realization of ferroelectric 
(ferromagnetic in the case of magnetic materials) phase, 
dipole (spin) glass phase and mixed ferroglass phases. 
Also, so-called paraglass (Griffiths) phase may occur in 



disordered dielectrics (see, e.giiS). 

On the other hand, both ordered and disordered di- 
electrics have intrinsic nonlinearities, consisting in, e.g., 
dielectric hysteresis. The "interaction" between these 
nonlinearities and random fields lead to their renormal- 
ization so that the distribution function of random fields 
will also include a nonlinear contribution of random 
fields. Such calculation had been carried out in a dif- 
ferent context in Refpii. 

The calculations in RefJi incorporate a self-consistent 
dependence of distribution function of random fields on 
the third order nonlinearity coefficient a 3 . It has been 
shown that when the nonlinear coefficient is sufficiently 
large and positive, the results are strongly different from 
those in linear caset In particular, for 0:3 > 0, the phase 
transition temperature exceeds its mean field asymptotic 
value, while for 0:3 < the results are qualitatively the 
same as in linear case. As we have found out earlier, 
this phenomenon is due to the "generation of order by 
disorder" (or more precisely, a specific positive feedback 
generated by positive nonlinearity) taking place at 0:3 > 
0. 



II. GENERAL FORMALISM 

A. The distribution function of random fields 

Here we briefly review the main facts about the shape 
of distribution function of random fields with respect to 
nonlinear effects, a more detailed discussion can be found, 
e.g. in RefsLii for disordered dielectrics and in RefsiS*^ 
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for disordered magnetic semiconductors. 

The distribution function of random field E can be 
represented in the form 



f{E) = (5(E-E{? i ))) 



(1) 



Here the bar denotes averaging over spatial configura- 
tions of random fields sources (e.g. electric dipoles, "re- 
sponsible" for emergence of ferroelectricity in the above 
compounds), (...) means the averaging over dipoles ori- 
entations, E(fi) is the internal electric field induced by 
electric dipoles and other sources in the observation point 
r, . In a disordered ferroelectric this field already contains 
the intrinsic nonlinearity and can be written in the form 
(see— for details) 



E 1 (r i ) =£ 1 {r i )+ ^ a m ]^[ £ 7j (r,;), 



(2) 



where 7=x,y,z so that £ 7 denotes simply the 7- 
component of vector £. Here £ is an internal electric 
field induced by electric dipoles, a m is a coefficient of 
nonlinearity of rn-th order. Note that the first term in 
Eq. can be generalized to account for other possi- 
ble sources of random fields such as point charges, elastic 
dipoles et o 5 ! 7 ! 11 . 

The calculation of the distribution function with 
respect to in the framework of statistical theory 
(se o 14 i 15 for details of this theory) for so-called disordered 
Ising model (when the dipole has only two admissible 
orientations) yields the following rigorous result for any 
electric field component E a = 

/oo / 00 \ 

h(E')5 [E-E'-Y, <*mE' m dE' (3) 
-°° \ ra=2 J 

where fi(E) is the distribution function, that takes into 
account only the first linear term in © (distribution 
function of the first order). 

More detailed version of Eq. © reads 

f(E)= [°° h(E')x 
x6 Ie - E' -jr a m E' m j dE', 

\ m=2 ) 

fi(E) = ±- r e iE *- nF Udt, 

27T J-oo 

■(( 



F(t) 



d 6 r 1 



-it£(r] 



(4) 

(5) 
(6) 



Here n and £(r) are the concentration and electric field of 
the dipoles. The Eq.© determines the function fi(E), 
calculated earlier ini for the case of two-orientable elec- 
tric dipoles. We note here, that in general case (e.g. ar- 
bitrary interaction between above dipoles, their arbitrary 
concentration etc, see Refi£ and references therein for dis- 
cussion) function fi(E) has non-Gaussian form. It can 
be shown that for our case of compound ferroelectrics it 
is sufficient to use its Gaussian asymptotics, which reads 



HE) 



1 



exp 



(E - E a Lj 2 
AnB 



(7) 



Here L — ^J- and Eq = Air(nd* 2 ) / 'eq are the order pa- 
rameter (number of coherently oriented impurity electric 
dipoles or dimensionless spontaneous polarization) and 
the mean value of random field of electric dipoles (in the 
energy units), d* = |d7(eo — 1) is the effective electric 
dipole moment, 7 and £0 are, respectively, Lorentz factor 
and static dielectric permittivity of the host lattice, n is 
the concentration of electric dipoles. Coefficient B deter- 
mines the width of the distribution function (dispersion 
of random fields) and depends on host lattice parameters 
like its correlation radius r c , see 7 - for details. 



B. The equation for long-range order parameter 

In our approach, an average value A of any physical 
quantity can be represented in the form 



A = 



f(E)A(E)dE, 



(8) 



where f(E) is determined by (@J and A(E) is the above 
quantity for single random field constituent, averaged 
over its internal degrees of freedom. In our case this is 
the average value over the orientations of single electric 
dipole. 

To calculate the T c for disordered ferroelectric com- 
pound we need to calculate first the long-range order pa- 
rameter L. In the spirit of Eq. (JHJ we obtain following 
self-consistent equation for this parameter 



fi{E) tanh 



E 



a m E m /k,T 



dE. 



(9) 
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Here we use the fact that for two-orientable dipoles (l z 
±1, l x = l y = 0) A{E) =< I >= tanh (J,), E = E z . 



The self-consistency of Eq. (jSJ is revealed by substitu- 
tion of Eq.lJSJ into it, which yields 
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exp 



(E - E Q L) 2 
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In the linear case (a m = 0) Ea. H10|) transforms into 
that derived in£. It is seen that order parameter is self- 
consistently expressed through itself and is the function 
of temperature, dipoles concentrations and nonlinearity 
coefficients. 

Further simplifications of Eq. (|1(J|> are possible on sym- 



J 



tanh 




dE. 



(10) 



metry grounds. Namely, for the lattice with the center of 
inversion in paraelectric phase, the order parameter has 
to be an odd function of electric field, i.e. m's in Ea. Hl(J|) 
are odd numbers. Conserving only the first nonlinear 
term in tanh argument, we obtain 
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exp 



(E - E LY 
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tanh [(E + a 3 E 3 ) /kT] dE. 



(11) 



r 



It is now instructive to consider the mean field limit of 
equation (jlip. This limit corresponds to the case of an 
ordered ferroelectric, where the distribution function of 
random fields degenerates into a delta function S(E — 
EqL). Formally, in our method this limit corresponds to 
7irj? — > oo. Substitution of this <5-function into Hll[) gives 
the desired mean field equation for order parameter 



J mf 



tanh 



^ cm/ 



{L m f + aoi m /) 



(12) 



where Eq = kT cm f, T cm f is transition temperature in a 
mean field approximation (see below), ao = ^Eq 80 that 
ao is dimensionless. We now use the Eq. (|ll|l to derive 
the equation for ferroelectric phase transition tempera- 
ture beyond the mean field approximation. 



C. Ferroelectric phase transition temperature 

The phase transition temperature is defined as a tem- 
perature when a nonzero order parameter L appears. In 
other words, to get an equation for T c from (|11|) . we 
should put in it L — ► 0. This is accomplished by noting 
that at small L 



exp 



Cg - EqL) 
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(13) 



Subsequent substitution of l|13|) into l|ll|) yields after 
some algebra 



A f (1 + 3a a; 2 )e-T ; 



o 



cosh 



2 X + QiQX^ 



-dx = 1. 



(14) 



Here wc introduced the following dimensionless variables 



A = 



Eq 
■kjiB 
E 



= ^/l5w;? 
kT c 



T, 



cmf 



(15) 



Equation (jl4|l is the main theoretical result of this work. 
It predicts the existence of the critical concentration of 
dipoles 

cr c ^ 15 ' 

such that for n < n cr , the long range order in the system 
would never be realized. Thus, the critical concentration 
is determined from the condition, that A = A cr at r c = 0. 
Taking the limit r c — > Oii in Eq. I|14|l , we obtain A cr = 1 , 
which justifies the choice of dimensionless parameter A. 
It is seen that the critical concentration does not depend 
on the coefficient of nonlinearity and it is completely the 
same as in linear case (se oV 1 ). 

Now we demonstrate that in a mean field approxima- 
tion r cm f = 1, i.e. that Eq — kT cm f. For that we notice 
that this approximation corresponds to A — * oo in Ijl4|l . 
which after some transformations^ gives T cm f = 1 . This 
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FIG. 1: Theoretical dependence r c (A 1 ). Numbers near 
curves denote magnitudes of non-linear coefficient, ao- 

value also does not depend on the coefficient of nonlin- 
earity. 

The plot of dependence r c (A _1 ) at different ao is shown 
on the Fig^ It is seen that at ao > the dependence 
t c (A -1 ) has maximum, while at ao < it does not. More- 
over, for negative nonlinearity there is sharp lowering 
(but not to zero) of r c at certain A -1 . This demonstrates 
that a negative feedback almost destroys long-range or- 
der especially at small dipoles concentrations. This be- 
havior shows that at ao > nonlinear effects produce 
positive feedback thus enhancing the long-range order in 
the system, while ao < the feedback is negative so 
that long-range order is inhibited but not completely de- 
stroyed even for large negative ao, see curve for ao = — 20 
on the FigQ] Also, at sufficiently large positive ao we can 
achieve very large enhancement of T c as compared to its 
mean field value (see curve for ao = 3 on the FigUJ. 
Since in our model we suppose that the value of T cm f 
is equivalent to T c for pure PbTiC>3 (the end member of 
compound), we conclude that large positive ao give sub- 
stantial increase of T c as compared to the end members 
of ferroelectric compound. 



III. COMPARISON WITH THE EXPERIMENT. 
DISCUSSION. 

In our model, the ferroelectric compounds are con- 
sidered as an ensemble of electric dipoles embedded in 
some virtual paraelectric host, its nature we will dis- 
cuss later. It can be supposed that in the considered 
system (BiSc03)i_ x (PbTi03) x electric dipoles are origi- 
nated from PbTiC>3, i.e. their number increase with x in- 
creasing. The nondipolar random field sources like point 
charges, elastic dipoles etc. are also present in such com- 
positions due to the mixed valency of Bi and the differ- 
ence in charges and ionic radii of Bi 3+ and Pb 2+ , Ti 4+ 
and Sc 3+ . These defects could easily be incorporated in 
the consideration (seeii for details), but we do not con- 



sider their contribution here to have the minimal number 
of adjustable parameters when fitting the theory with ex- 
periment. 

For quantitative description of experiment in above fer- 
roelectric compounds in the framework of our random 
field theory (in linear or nonlinear approximation) we 
should have precise information about the concentrations 
of electric dipoles and other random fields sources as well 
as the parameters £q, r c etc. Unfortunately the avail- 
able data are strongly restricted for above compounds. 
That's why we recalculated our dimcnsionless parameters 
(|15|l from the best fit to the experiment^. Namely, from 
the position and " amplitude" of maximum of experimen- 
tal curve T c (x) we determine our nonlinearity coefficient 
ao = 0.81 and coefficient of recalculation of parameter A 
into x. Latter coefficient gives us value T cm f = 490°C 
and critical content of PbTiC>3 x cr — 32%. The result of 
such fitting is shown on the FigEl 
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FIG. 2: Comparison of our theory (line) with experiment- 
(squares). Inset shows the fit in more details. 

It is seen in Figure 2 that there is pretty good coinci- 
dence between theory and experiment. This coincidence 
shows, that physical mechanism of enhancement of fer- 
roelectric phase transition temperature beyond the end 
members in the ferroelectric compounds is the enhance- 
ment of initial (i.e. that of end members) long range 
ferroelectric order in them caused by positive feedback 
generated by the nonlinear effects with positive coeffi- 
cient. The obtained value T cm t = 490° C is very close to 
transition temperature of PbTiC>3. This speaks in favour 
of statement that PbTiC>3 paraelectric phase can be con- 
sidered as a host lattice. 

As we have mentioned above, the origin of nonlinear ef- 
fects in compound ferroelectrics is quite naturally related 
to its intrinsic nonlinearities, reflected, for instance, in 
their hysteresis loop. The physical origin of these nonlin- 
earities may be the nonlinear coupling of the ions in the 
unit cell as well as the clusterization of dipoles in diluted 
ferroelectric compounds. Each cluster has its own meso- 
scopic dipole moment (polarization of cluster) and such 
clusters can interact between each other. If interaction 
between such clusters is of positive sign, which may be 
manifested in the almost rectangular shape of hysteresis 
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loop, we have the positive feedback with enhancement of 
initial (i.e. that of end member) long-range order with 
positive nonlinear coefficient. If the interaction is of neg- 
ative sign, we have negative feedback and inhibition of 
ferroelectric order. But such negative feedback cannot 
destroy ferroelectricity completely, probably except for 
the case of extremely small dipoles concentration. So, 
to achieve high positive values of nonlinear coefficients 
(which are necessary to strongly enhance T c ), we should 
have (or prepare) the end members of ferroelectric com- 
pound with "as much as possible rectangular" hystere- 
sis loop. The contribution of nonlincaritics to dielectric 
response and hysteresis loop in mixed perovskites PZT 
with impurities was discussed in Refii£. The essential 
role of nonlinearity effect in explanation of phase dia- 
gram peculiarities in mixed ferroelectrics PbZri-^Ti^Os, 



BaZr^Tii-^Oa and mixed systems of ferroelectric relax- 
ors was shown recently in Refi^ and Ref»2S respectively. 

To make more precise prediction of which components 
in ferroelectric compound to use to increase T c beyond 
end members value, the measurements of correlation ra- 
dius, Lorentz factor and ions shifts in both end members 
and entire above ferroelectric compounds arc highly de- 
sirable. 

It should be finally noted, that the theory outlined 
here permits one to calculate not only T c for above sub- 
stances but virtually any equilibrium and/or nonequilib- 
rium thermodynamic characteristics such as static and 
dynamic dielectric susceptibility, specific heat etc as a 
functions of PbTiC>3 content, temperature, electric field 
and other external parameters. 
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